Search for: All records

Variable viscosity in Earth’s mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 $\pm$0.25 (mean $\pm$SD). A relatively low yield stress of 151 $\pm$19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes. 

Abstract Growth of the Andes has been attributed to Cenozoic subduction. Although climatic and tectonic processes have been proposed to be first-order mechanisms, their interaction and respective contributions remain largely unclear. Here, we apply three-dimensional, fully-dynamic subduction models to investigate the effect of trench-axial sediment transport and subduction on Andean growth, a mechanism that involves both climatic and tectonic processes. We find that the thickness of trench-fill sediments, a proxy of plate coupling (with less sediments causing stronger coupling), exerts an important influence on the pattern of crustal shortening along the Andes. The southward migrating Juan Fernandez Ridge acts as a barrier to the northward flowing trench sediments, thus expanding the zone of plate coupling southward through time. Consequently, the predicted history of Andean shortening is consistent with observations. Southward expanding crustal shortening matches the kinematic history of inferred compression. These results demonstrate the importance of climate-tectonic interaction on mountain building. 

Abstract The dip angles of slabs are among the clearest characteristics of subduction zones, but the factors that control them remain obscure. Here, slab dip angles and subduction parameters, including subduction duration, the nature of the overriding plate, slab age, and convergence rate, are determined for 153 transects along subduction zones for the present day. We present a comprehensive tabulation of subduction duration based on isotopic ages of arc initiation and stratigraphic, structural, plate tectonic and seismic indicators of subduction initiation. We present two ages for subduction zones, a long‐term age and a reinitiation age. Using cross correlation and multivariate regression, we find that (1) subduction duration is the primary parameter controlling slab dips with slabs tending to have shallower dips at subduction zones that have been in existence longer; (2) the long‐term age of subduction duration better explains variation of shallow dip than reinitiation age; (3) overriding plate nature could influence shallow dip angle, where slabs below continents tend to have shallower dips; (4) slab age contributes to slab dip, with younger slabs having steeper shallow dips; and (5) the relations between slab dip and subduction parameters are depth dependent, where the ability of subduction duration and overriding plate nature to explain observed variation decreases with depth. The analysis emphasizes the importance of subduction history and the long‐term regional state of a subduction zone in determining slab dip and is consistent with mechanical models of subduction. 

Abstract The extensive fast seismic anomalies in the mantle transition zone beneath East Asia are often interpreted as stagnant Pacific slabs, and a reason for the widespread tectonics since the Mesozoic. Previous hypotheses for their formation mostly emphasize vertical resistances to slab penetration or trench retreat. In this study, we investigate the origin of these stagnant slabs using global‐scale thermal‐chemical models with data‐assimilation. We find that subduction of the Izanagi‐Pacific mid‐ocean ridge marked the transition of mantle flow beneath western Pacific from being surface‐driven Couette‐type flow to pressure‐driven Poiseuille‐type flow, a result previously unrealized. This Cenozoic westward mantle wind driven by the pressure gradient independently explains seismic anisotropy in the region. We conclude that the mantle wind is the dominant mechanism for the formation of stagnant slabs by advecting them westward while the pressure gradient holds them in the transition zone.